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ABSTRACT 

We generalize Starobinskii's stochastic technique to the theory of a massless, 
minimally coupled scalar interacting with a massless fermion in a locally de 
Sitter geometry. The scalar is an "active" field that can engender infrared 
logarithms. The fermion is a "passive" field that cannot cause infrared loga- 
rithms but which can carry them, and which can also induce new interactions 
between the active fields. The procedure for dealing with passive fields is to 
integrate them out, then stochastically simplify the resulting effective action 
following Starobinskii. Because Yukawa theory is quadratic in the fermion 
this can be done explicitly using the classic solution of Candelas and Raine. 
We check the resulting stochastic formulation against an explicit two loop 
computation. We also derive a nonperturbative, leading log result for the 
stress tensor. Because the scalar effective potential induced by fermions is 
unbounded below, back-reaction from this model might dynamically cancel 
an arbitrarily large cosmological constant. 
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1 Introduction 



Massless, minimally coupled scalars and gravitons are unique in achieving 
zero mass without classical conformal invariance. This means that inflation 
rips their virtual quanta out of the vacuum, which greatly strengthens the 
quantum loop effects they mediate [I]. In the expectation values of famil- 
iar operators these enhanced quantum effects typically manifest as infrared 
logarithms. A simple example is provided by the stress tensor of a massless, 
minimally coupled scalar with a quartic self-interaction, 

C = —-^d^j,(pdyipg^" — ^f^^y/^ + counterterms . (1) 

When the expectation value of the stress tensor of this theory is computed 
in de Sitter background, 

ds^ = -dt^ + a^{t)dx ■ dx with a{t) = e^* , (2) 

and renormalized so as to make quantum effects vanish at t = 0, the results 
for the quantum- induced energy density and pressure are [21 13] , 

Pit) = (3) 

Infrared logarithms are the factors of ln(a) = Ht. They arise from the fact 
that inflationary particle production drives the free scalar field strength away 
from zero 14], El E] 

if'^ix) ^ij^ = ln(a) -|- Divergent Constant . (5) 

This increases the vacuum energy contributed by the quartic potential and 
the result is evident in 

Infrared logarithms arise in the one particle irreducible (IPI) functions 
of this theory |7] . They occur as well in massless, minimally coupled scalar 
quantum electrodynamics (SQED) [H [9l [TOl [11] and in massless Yukawa the- 
ory [T2l|13]. The IPI functions of pure gravity on de Sitter background show 
infrared logarithms [Tll[T5l[IS]. A recent all orders analysis of scalar-driven 
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inflation was unable to exclude the possibility that they might even contami- 
nate loop corrections to the power spectra of cosmological perturbations [T7] ! 
And infrared logarithms have been discovered in the IPI functions of Gravity 
+ Dirac 

Infrared logarithms are fascinating because they introduce a secular el- 
ement into the usual, static expansion in the loop counting parameter. No 
matter how small the coupling constant A is in the continued growth of 

the inflationary scale factor must eventually overwhelm it. When this hap- 
pens perturbation theory typically breaks down. For example, the general 
form of the induced energy density ([3]) is, 

p{t) = H'p^X'-'I^C',ln''-\a) + Ciln''-\a) + ... + C',,_,\n\a)'^ . (6) 

The A^Cq ln^^~^(a) terms are the leading logarithms at i loop order; the re- 
maining terms are suhdominant logarithms. Assuming that the numerical 
coefficients are of order one, we see that the leading infrared logarithms 
all become order one at ln(a) ~ l/\fX. At this time the highest subdominant 
logarithm terms are still perturbatively small (~ \/A), so it seems reasonable 
to attempt to follow the nonperturbative evolution by resuming the series of 
leading infrared logarithms, 

oo 

Plead = i^'E^o(Aln'(a)) ~ . (7) 

This is known as the leading logarithm approximation. 

Starobinskii has long maintained that his stochastic field equations repro- 
duce the leading logarithm approximation [2D] . With Yokoyama he exploited 
this conjecture to explicitly solve for the nonperturbative, late time limit of 
any model of the form |21j . 

C = -U,^d,^g^''^-V{^)^ , (8) 

assuming only that the potential V{ip) is bounded below. When the potential 
is unbounded below the conjecture still gives the leading infrared logarithms 
at each order, however, the theory fails to approach a static limit. 

Starobinskii's conjecture has recently been proved to all orders (221 [23] . 
The field equations are first rewritten in Yang-Feldman form [21], then the 



2 



free field mode expansion is truncated at horizon crossing, and tlie free field 
mode functions are replaced with their leading long wave length forms. This 
procedure converts the original quantum field into a commuting random vari- 
able, but it preserves the leading infrared logarithms. Although it was not 
evident at the time, the reason the field can be infrared truncated is that 
every pair of fields in a simple potential model of the form ([HD is capable 
of inducing an infrared logarithm, and the leading log term derives from 
requiring them to do so. Because only the infrared part of the field is respon- 
sible for infrared logarithms, we can truncate and take the long wave length 
limit of the mode functions, For example, the infrared 

logarithm in ^ is, 

A field which can generate infrared logarithms is called active. Scalar 
potential models of the form possess only active fields. However, more 
general theories can possess fields which are not themselves capable of en- 
gendering an infrared logarithm. We call these passive fields. A example of 
such a model is SQED, 

C = -^-F^^F.^g^Pg^'^^ - (d,-zeA,y (d,+zeA,)^g>^''^ . (10) 

In this model the charged scalar is active whereas the photon is passive. 




Fig. 1: Two loop contribution to {Q\F^^{x)Fp„{x)\Q) . 

Although passive fields cannot cause infrared logarithms, they can prop- 
agate their effects. That is, an expectation value of passive fields can acquire 

^It is also necessary to take the first nonzero term in the long wavelength limit for 
the retarded Green's function of the Yang-Feldman equation. This turns out to require a 
higher order term in u(t, k). 
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an infrared logarithm from a loop correction involving an active field. For 
example, the diagram in Fig. 1 gives a contribution to {fl\F^^{x)Fpi^{x)\fl) 
which acquires an infrared logarithm through the scalar loop at the bottom. 

Passive fields can also induce interactions between active fields. For ex- 
ample, the photon loop in Fig. 2 induces an effective [ip*ifY interaction in 
SQED, 




Fig. 2: Effective {(p*(p)'^ coupling in SQED. 

SQED — and also gravity — feature another complication in which 
derivatives of active fields can induce interactions between undifferentiated 
active fields. For example, the 3-point interaction of SQED, 

te{ip*d^^ - d,ip*ip)A,g^''^ , (11) 

can induce an effective (p*(p coupling through the diagram of Fig. 3, 




(f*d(f dip* if 
Fig. 3: Effective (p*ip coupling in SQED. 

This is part of the full IPI 2-point function which has recently been computed 
at one loop order [25] . 

In generalizing Starobinskii's technique to theories which include passive 
fields, and/or differentiated active fields, it is crucial to realize that the ultra- 
violet parts of passive fields and differentiated active fields contribute on an 
equal footing with the infrared parts in propagating infrared logarithms and 
in mediating interactions between undifferentiated active fields. So one can- 
not infrared truncate the passive fields, or even differentiated active fields. 
Instead the correct procedure is: 
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1. Integrate out the passive fields and renormalize the resulting effective 
action. 

2. Integrate out the differentiated active fields and renormalize the result- 
ing effective action. Note that this can always be done because the 
original action is at most quadratic in derivatives. 

3. Infrared truncate and stochastically simplify the effective action of ac- 
tive fields. 

One might suspect that the third step is not possible owing to the nonlocal- 
ity of the effective action. However, this nonlocality must be mediated by 
differential operators which, precisely because they derive from inactive fields, 
must contain positive powers of the scale factor whose rapid time dependence 
weights the integral overwhelmingly at its upper limit and totally dominates 
the logarithms which might derive from the active fields. A typical example 
is the integral, 

[' dt'a'^\n{a') = [' dt'e^""'' Ht' = — |ln(a') - -| 
Jo Jo 2H ^ 2 J 

For In (a) ^ 1 it is as though we simply divide the integrand by 2H and 
evaluate it at the upper limit! Hence the hopelessly complicated "effective 
action" degenerates, in the leading log approximation, to a very tractable 
"effective potential," and the resulting local theory assumes the form ([8]) 
already solved by Starobinskii [201 [2T]. 

Yukawa theory is especially simple because it possesses no differentiated 
active fields, and because it is free of the subtle gauge fixing problems of 
SQED [25j and gravity f2E{ \TE\ . In section 2 we review the full apparatus of 
perturbation theory for massless Yukawa on a locally de Sitter background. 
In section 3 we integrate out the fermion and renormalize the effective poten- 
tial. A curious and possibly significant property of Yukawa theory is that its 
effective potential is unbounded below, a fact that survives in the flat space 
limit and has even earned a place in Standard Model parameter estimation. 
We check the stochastic formalism against an explicit two loop vacuum ex- 
pectation value in section 4. In section 5 we employ the stochastic formalism 
to obtain the leading log approximation for the Yukawa stress tensor. Be- 
cause the effective potential depends upon the inflationary Hubble constant, 
the induced vacuum energy of the stress tensor does not quite agree with it. 



^_La2in(a). (12) 
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In fact the latter is initially positive whereas the former is always negative. 
However, their asymptotic large field behaviors are identical, so there seems 
no avoiding the conclusion that gravitational back-reaction in this model 
must eventually halt inflation, although not in an acceptable fashion. Our 
conclusions comprise section 6. An appendix presents the less interesting, 
technical analysis behind the perturbative calculation described in section 4. 



2 Massless Yukawa in de Sitter 

The coupling of gravity to particles with half integer spin is usually accom- 
plished by shifting the fundamental gravitational field variable from the met- 
ric g^iu{x) to the vierbein e^m(x), although there are other approaches [27] . 
Greek letters stand for coordinate indices, Latin letters denote Lorentz in- 
dices, and both sorts take values in the set {0, 1, 2, ... , {D—\)}. One recovers 
the metric by contracting two vierbeins into the Lorentz metric rf'^, 

Qiiuix) = e^,b{x)e^c{x)rf''' ■ (13) 

The coordinate index is raised and lowered with the metric {e^ b = g^'^Gub), 
while the Lorentz index is raised and lowered with the Lorentz metric (e^ ^ = 
rf'^e^c)- We employ the usual metric-compatible and vierbein compatible 
connections, 

gpa;t, = ^ = ]^g'"' {ga^,,u + gua,t, - gpu,a) , (14) 

e/36;M = =^ Af,cd = e''c(e^d,ti - '^^uepd) ■ (15) 

Fermions also require gamma matrices, 7^^. We assume their spinor in- 
dices run over i = 1, . . . , 4 in any dimension. The anti-commutation relations 
are, 

{7^ 7'=} ^ (7^^ + 7^7") = -2rj'^I . (16) 
The Dirac Lorentz representation matrices are. 



jbc — ^ 

^ -4 



7 ,7 



(7V - 7^7') • (17) 



They can be combined with the spin connection (llSp to form the Dirac co- 
variant derivative operator, 

V^ = d^+ '-A^cdJ"" . (18) 



6 



In a general vierbein background the bare Lagrangian of massless, mini- 
mally coupled, Yukawa scalars with massless fermions is, 

+ii,e\^''Vp^^ - hv'^i^^ . (19) 

The bare fields in this expression are ^p{x), and il)i{x). The symbols 

^0, Ao and /o stand for, respectively, the bare conformal coupling, the bare 
4-point coupling constant and the bare Yukawa coupling constant. Neither 
the scalar nor the fermion requires a mass term because we desire the special 
model with zero renormalized masses, and no mass counterterms are required 
because mass is multiplicatively renormalized in dimensional regularization. 

Renormalization is begun by expressing the bare fields in terms of the 
renormalized ones, 

(y9 = ^fZ i^y. , l\} = \J^2'^r ciud ijj = \fZ2lpr: ■ (20) 

Substituting fl20l) into the bare Lagrangian f|T9|) gives, 

+tZ2^y,'y''Vp^,,/^ - y/ZZ2fo^,'^r^rV^ ■ (21) 

We next enforce the conditions that the renormalized scalar should have 
neither a conformal coupling nor a 4-point coupling, 

Z^o = + 6^0 , Z^Xo = + 6X0 and v^Za/o = / + 5/ • (22) 

Of course the model does require conformal and 4-point counterterms! We 
also define the field strengths as usual, 

Z = 1 + 6Z and Z2 = 1 + . (23) 

The structure of renormalized perturbation theory is complete when we ex- 
press the Lagrangian in terms of primitive interactions and counterterms, 

- yV^^^V^ - - ^M^r^rV^ • (24) 
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The preceding analysis has so far been for a general geometry. Of course 
we are interested in the special case of inflation, for which de Sitter is an 
excellent paradigm. We work on the open submanifold of D dimensional 
de Sitter space in the conformal coordinate system for which the invariant 
element is, 

ds^ = a?ir]){—d-rf + dx-dx) and air]) = . (25) 

^ ^ Hr] 

Of course this makes the metric g^^ = a^{vi)rj^^. We also employ the ghost- 
free, Lorentz symmetric gauge [28| which determines the vierbein, 

e^b = ebf, =^ Cf^b = ar]f,b ■ (26) 

At this stage there is no more point in distinguishing Lorentz indices from 
coordinate ones. With these conventions the covariant derivative operator 
takes the simple form, 

V,-^d, + ^Ha[^',j,] . (27) 

The special case of its contraction into e^^7^ is even simpler, 

^'e^^V.^i;^ a-(^)^a(V) ^ a-^'^h^d.a^'^^ . (28) 

The scalar and fermion propagators can be largely expressed in terms of 
the following function of the invariant length i{x; x') between and x''^, 

yix;x') = 4:sm'^(^^He{x;x')) = aa'H^Ax'^{x;x') , (29) 
= aaH^{\\x-x'f -{\7]-7]'\-t6y) . (30) 

The most singular term for each case involves the propagator for a massless, 
conformally coupled scalar, 

^Acf(x;x') = ^r(^-l)(-f'' . (31) 



(An) 



2 



V 2 



It has long been known that there is no de Sitter invariant solution for the 
propagator of a massless, minimally coupled scalar [29j. If one elects to 
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break de Sitter invariance while preserving homogeneity and isotropy — this 
is known as the "E(3)" vacuum [20] — the minimal solution is [H [3], 

iA{x; x') = iAcf{x] x') 



+ 



H^~' V{D-l] 



D 



I 2 



(47r 

ttD-2 oo 



[D-4:T{D-1) 



4xf-2 



TTcoty—Dj + In(aa') 



'1 T{n+D-1) fV 



(47r)-^iln T{n + 



r(n + f + l)^?/^n-f+2| 



■(f)' 



.(32) 



n-f + 2 r(n + 2) 

This expression may seem daunting but it is actually simple to use because 
the infinite sum vanishes in D = 4, and the terms of this sum goes like higher 
and higher powers of y{x;x'). Hence the infinite sum can only contribute 
when multiplied by a divergent term, and even then only the first few terms 
can contribute. 

It is useful, in the stochastic analysis to follow, for us to consider the 
fermion propagator in the presence of an arbitrary, potentially spacetime 
dependent, scalar field ip{x), 



iS[f^p]{x;x') = (x 



X 



(33) 



iS[m]{x; x) 



+ \l —ml 



Of course this can only be evaluated for a handful of field configurations. 
The case of a constant, fip{x) = m, was solved by Candelas and Raine [31j . 

^ - 1 1 ( at'U> - ^ 

D 

' D \ 

~2 



(47r)T V9 A ^ 



'aa' 



r(f-i+zf)r( 



X ■ 



2 * 



r( 



r(f-i)r(,. 



i-zf)r(f+zf; 



-2i^i(--i+v.> 



m D ,m D 



)( 



r(f-i)r(f) 



.m D 



2Fi{--l-t—,-+t—---l--^ 



,m D 



4 



(34) 



The propagator we use in perturbative calculations is the massless limit of 
this one, 



S 



--^r(--i a^^ 

(47r)- ^ 2 ^ 
ailp —j= iAcf (x; x') 



'aa 



D D y 
1 — • — • 1-- 

' 2' 2' 4 



[aa 



aa' 

-1 

2 ' X t( 



X 



r(f-i) 



47rT • 



(35) 
(36) 
(37) 



The final expression is just a conformal rescaling of tlie propagator for a 
massless fermion in fiat space, as it sliould be in conformal coordinates. 

It is useful to recast the Candelas-Raine solution flM|) using the transfor- 
mation formula of hypergeometric functions (see expression 9.131.2 in [32j). 



iS[m]{x;x') = ,, r(—-l)(ai'P 



' — ml 

a' J 



(47r)T ^2 ^ ^ ' y/aa' V a 

r(i-f)r(f-i+zf)r(f-zf) ^jD_ 

2 ^~^'H' 2 2 ' 4 

D //-7° 



r(f-i)r(^f)r(i-zf 



/ly m D m D y 
FA i—;—;- 



+ 



.m .m ^ 
I — ,1 — ;2 — - , 
H' H' 2 '4' 



+ 



r(i-f)r(f-i-zf)r(f+^g 
r(f-i)rHf)r(i+^ 



m D m D y 



4) '^'i'^'H^-'H'^^-r-AJ V 2 



if' 2 H' 2' A' 



In dimensional regularization all D-dependent powers of y{x\ x') vanish at 
coincidence. Hence we obtainj^ 



V -cr u /^ mii^-'r(f+zf)r(f-if)r(i-f) ^ 



lim Vf^iS[m]{x; x') 
lim iS[m]{x] x')V ^ 



r(i+zf)r(i-^f) 
mii^-2r(f+zf)r(f-^f)r(i-f) 



(39) 



(4vr)f r(i+zf)r(i-^f) 
mg^-^r(f +zg)r(f-zg)r(i- 
(4vr)f r(i+zf)r(i-^f) 



X -^rna-i^ , (40) 



X 



— ma7^. (41) 



3 Stochastic Effective Action 

Integrating out the fermions gives rise to a scalar effective action, 



e'^[^''^''^l = e'^=[^l det[v^(^?* - hv 



(42) 



^Note that the back-acting covariant derivative is V ^ = 9^ — |i/a[7",7^]. 



10 



Here Ss[^p] stands for the action associated with the purely scalar part of the 
bare Lagrangian, 



= -U.^d^vg^'^^g - %v'R^ - ^v'V^ . (43) 
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The exact effective scalar field equation is, 



-Tr 



5 



-9{iV -fo^) S^{x) 
An 



-9[^ 



d.^^g^-'dM^)) - iM^)R^g - -^ip\x)^g 



+/oTr iS[foip\{x;x) 



-9 



,(44) 



(45) 



Were we trying to solve the full quantum field theory, equation fH^ would 
be a dead end because we lack an explicit expression for the coincidence limit 
of the fermion propagator in the presence of a general ^p{x). However, we are 
focused instead on the leading infrared logarithms, and this fact permits a 
crucial simplification: we can evaluate the fermion propagator as if the scalar 
were constant! Making use of (1391) we infer a local field equation which agrees 
exactly with the full theory in the leading log approximation, 



6T\ 



6(p{x) 



-gg^^d,if]-^o^R, 



X 



-9 



-9 



4/0^/7^-2 r(f+.^; 



r(f-.if)r(i- 



' 2 > 



(4vr: 



r(i+.if )r(i-.if 



■^^g. (46) 



Of course the final term in fj46|) is also related to the effective potential (it is 
—Vl^{ip)^—g^ and has appeared many times before in this guise [31| [33 l [34 



The factor of r(l — y) in equation fB6l) is the only divergence we shall 
see in the stochastic formalism. It can be removed using the parameters of 
the scalar potential, and Aq. In particular, the stochastic formalism does 
not require either field strength renormalization or renormalization of the 
Yukawa coupling. 



5Z = = (5/ = 



^0 = 5^ and Ao = 5A . 



(47) 



11 



To renormalize fl46l) we set D = 4— e and make use of the following expansions, 

(48) 



r(l-y)=-^+7-l + 0(e), 



H 



H 



+ 0{e^) . (49) 



The symbol "'?/'(z)" in this last expression of course stands for the psi function 
rather than the fermi field (see section 8.36 of [32j), 



d 



^{l+z) ^ - ln(r(l + ;.)) = -7 + E(-irC( 



dz 

Note that fj49l) is both real and even in (p, 

r(f+^4^)r(# 
r(i+ 



n]z 



n—l 



(50) 



n=2 



- E(-ir [C(2n-1)-C(2n + 1)] (^) "e + 0{^) . (51) 



From these expansions it is apparent that we can renormalize so as to keep 
the scalar massless and free to order 



2 7TD-4 



5e 
5X 



r 1-4^ 



^ + :e?-^(l-7) 



(47r)f D(D-l) 247r2 



24/^17 



4 ijD-i 



(4vr: 



-r 1 



+ 



3/\ 



C(3)-7 



(52) 
(53) 



Substituting in fH6|) and taking the limit D A gives the following effective 
equation of motion, 



2^2 E(-l)1C(2^-l)-C(2^ + l)](^) " v^, (54) 



X 



T?C(3) 

27r^ 
1 

2' 



if 



27r2 
if 



(1)^35. ,55, 
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It is well to digress at this point to establish an important correspondence 
limit that bears upon the validity and physical interpretation of our renor- 
malization condition fl52l) . Duffy and Woodard computed the one loop scalar 
self-mass-squared in this theory and used it to solve the effective scalar field 
equation |35]. They found (equations (38) and (77) in [35]) that the following 
choice for the conformal counterterm results in there being no significant late 
time corrections to the scalar mode functions at one loop order, 



Si 



DW 



■D 
2 ' 



(47r)f 2(Z}-l)(D-3) 
(47r)T ^ 3e^ IGtt^ 



+ 



327r2 

1 
1 



Of 



(56) 



(57) 



Up to irrelevant terms of order e, this is precisely the same renormalization 
fl52|) as we have used! Therefore our stochastic renormalization conventions 
agree with the full theory in the regime of significant late time effects, just 
as they should. 

Up to some different renormalization conventions, our expression for the 
scalar effective potential agrees with equation (40) of Candelas and Raine 

ED, 




1)-C(2n + 1) 



(58) 



+2 J dx{x+x'^)\ij{l+ix)+ilj{l-ix)\j. (59) 

We have now reduced the theory to a completely finite, scalar model of 
the form already solved by Starobinskii [201 [21]. An interesting and possibly 
significant result of applying his technique is that this model fails to approach 
a static limit at late times! This is obvious once one recognizes that the 
potential V{ip) is unbounded from below. 

One might expect that V{(p) is negative because a nonzero scalar (of either 
sign) drives the fermion mass positive [121 [IS] 5 which must lower the vacuum 
energy. The absence of a lower bound is most easily proved by making use 
of the asymptotic expansion Stirling's formula implies for the psi function. 



m = Hz) - ^ - 



(60) 
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to evaluate the strong field limit of (15^ . 



Of course the large ip regime is also the small H regime, at which point 
we make contact with equation (6.10) in the classic paper by Coleman and 
Weinberg [3B]. The negative potential has long been recognized to render 
pure Yukawa unstable in fiat space [37j. A constraint on the Higgs mass can 
be derived in the Standard Model from the need to avoid this instability for 
the large Yukawa coupling of the top quark [38]. 

Although we have just seen that the instability is present in fiat space, 
inflation does have a role to play. A flat space scalar would simply roll 
down an unbounded, negative potential. However, the Hubble friction of 
expansion retards the scalar's downward progress. The potential's curvature 
is slight, for small /, so the scalar's evolution is for a very long time driven 
by the pressure of inflationary particle production. Only when the scalar's 
magnitude approaches the nonperturbatively large scale of 9? ~ H/ f does 
the unbounded potential begin to dominate. 



4 Reality Check 




Fig. 4: Lowest order contribution to {VL\T[(p^{x)il)^{x)il)r{x)]\Vl) . 



In this section we will test the stochastic formalism by comparing its pre- 
diction with an explicit two loop evaluation of the coincident vertex function 
(see Fig. 4), 



T 



iPr{x)ijj.{x)iJr{x) 



(62) 
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Because there is no field strength renormalization in the stochastic formalism, 
we can ignore the distinction between renormalized and unrenormalized fields 
in working out the stochastic prediction for fl62l) . Integrating the fermions out 
exactly gives the trace of the field-dependent fermion propagator, which we 
can again evaluate for constant field configurations in the leading logarithm 
approximation. 



iS[v,y^M X ip{x)i;{x)ip{x) = ^(x)Tr 



-iS[f(p]{x;x) 



H 



D-2 



X 



{4n 



r 1 



AfH 



D-2 



-r 1 



D 



D 



r(i+zf) 



0{f 



(47r)f ^ 2/ V2 
Hence our prediction for the order / infrared logarithm is divergent. 



n 



T 



(p,{x)lp^{x)lljr{x) 



n 



+ 



+ finite + 0(r 



(63) 
(64) 
(65) 

(66) 



»7r' e 

Note that the divergence arises from integrating out coincident fermion fields. 
It has nothing to do with the stochastic formalism per se, except for being 
the correct leading log result for the expectation value we have chosen to 
compute. 

That was easy! A measure of the power of the stochastic formalism is 
that we could just as simply have obtained the leading log result at order 

or higher. We turn now to the much more difficult task of perturbatively 
computing the full order / result for comparisonjl In models for which the 
"in" and "out" vacua differ, the in-out matrix elements computed with the 
usual Feynman rules are not true expectation values. To obtain an expecta- 
tion value such as (l62l) one must employ the Schwinger-Keldysh formalism 
[391 Hni im SSI HSl HH us] . For a recent review of the position-space formalism 
see |l6]. Here we simply summarize the modified Feynman rules: 

1. Each line has a polarity which can be either "+" or "— ". 

2. Vertices, including counterterms, are either all + or all — . 

3. A + vertex is the familiar one of the in-out formalism, whereas a — 
vertex is its complex conjugate. 



^This computation was done in collaboration with P. M. Ho. 
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4. External lines from time-ordered operators are +, whereas external 
lines from anti-time-ordered operators are — . 

5. Propagators can be ++, H — , — h or . For our theory these are 

all obtained from the Feynman propagators (1321) and (I57|) by replac- 
ing the conformal coordinate interval Ax^(x; x') with the appropriately 
polarized interval, 

i\v-v'\-^5)\ (67) 

{r]-r]'+t6)^ , (68) 
(^_^'_,5)2 ^ (69) 

{\r]-r]'\+i6f . (70) 



Ax^_^_{x; x') = \\x — x' 

Ax'^_{x;x') = \\x — x 

Axt^{x;x') = \\x — x 

Ax'^__{x;x') = \\x — x 



Because the operators in fl^ are all time-ordered, the associated external 
lines have -|- polarity. At lowest order they can connect to either a -|- or a — 
vertex. Hence the Schwinger-Keldysh result for (162|) is. 



n 



T 



iPr{x)lp^{x)lprix) 



+ifid''x'a''' 



n 

iA++{x] x')i 
-iA+_{x; x')i 



iSj 
iSj 



jSi 
jSi 



(X :x) 



0{f).{n) 



We will ignore the distinction between -|--|- and H — until it becomes signifi- 
cant. In view of (1371) the spinor trace gives. 



iSj 



jSi 



,x :x 



lijYjAx^Axy 



47r^ (aaO^-iAx2^ 



r^(f ) (ggQ^-^ 



(72) 



And the lowest order contribution to the VEV takes the form. 



V9r(a;)^/',(x)V'r(a;' 



-if 



d^x'a'^^^^^^^'^"^ ^A+_(a;;x') 



Axl^- 



+ Oif). (73) 



Here is where the analysis becomes tedious. We have accordingly moved 
the technical details to the Appendix, and will here quote only the diver- 
gent parts of the result. The superficially most singular contribution to fl73l) 
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derives from the first line of the scalar propagator We call it A, 

= finite . (75) 
The next most singular contributions come from the second line of fl32|) . 

_/i/V-^^(i-fi)(i-T)r^(2-f)r 1 Ha) 



{--^-^} + fimte, (77) 



x-.cot(^Z.)/.Va'{^-^), (78) 
-^fcot@r(l-.){l + M^} + fimte, (79) 

fH^-'fx-'.^ . r 3 ln(a). ^ . 

Only the n = 1 term from the third line makes a nonzero contribution in the 
limit of D = 4, 



CAx) = 



2 



D-6 W ^ i ^ ^ \A,t3^-8 Axl^-^j ^ ^ 
fH'f^-'' (l-i)(l-i)(l-fr(2-f) r 3 . 

^^^^^^ (l-|6)(l + e) 1^1 (83) 



^'('^)- 2^vrfV-i r(f + l) 4 lAx2^-4 Axi^-^r' (^4) 
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The = 2 term goes like Ax^, which is enough to make denominator inte- 
grable in D = 4. We can therefore take D = A for these terms, at which 
point one sees that they vanish. The same argument apphes to all terms 
with n>2. 

The sum of ([75]), ([77D, ([79]), ([HI]), ([83]) and ([85]) gives, 



( divergent constant ) + — — V unite , (86) 



in perfect agreement with the stochastic prediction ([66]) ! It is worth noting 
that the only the S-terms contribute divergent infrared logarithms, and those 
from B\ and cancel. So the result seems to derive entirely from i?3, which 
itself originated from the explicit factor of In(aa') on the second line of the 
scalar propagator ([32]) . 

Although it is not really necessary for our purpose of checking the stochas- 
tic formalism, we remark that the reason the coincident vertex diverges is 
that ordinary renormalization does not generally suffice for composite oper- 
ators such as ([B2]) . To renormalize local composite operators one must allow 
them to mix with operators of the same or lower dimension. There are three 
local, dimension four operators that can mix with >^^{^x^'i\} ^[x^'>\}^[x\ 

R^l , d^^,d,^,g^'' and i),e\-l'V^^\), . (87) 

The expectation value of the last term vanishes in dimensional regularization. 
The expectation value's of the first two are. 



T 



{a\T-\d,,p,{x)d,^,(x)g<-'\\a) = y^^T-i rp"^.' ■) (89) 

We can choose the coefficient of R^\[x) to completely cancel the second term 
of ( !86l) . Then we can choose the coefficient of dp,^^{x)dyif^{x)g^^ to cancel 
whatever constant terms remain. 
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5 Stochastic Stress Tensor 



To understand how this model sources gravitational back-reaction we must 
study the Yukawa stress tensor, 

1 6S 



T = 



-9 



(90) 



- ^^V • (91) 



Integrating out the fermions converts the fermionic terms to a purely scalar 
expression that we can evaluate for constant field configurations using (HOj) 
and (HID, 



la 

lim Tr 

2 x'^x 



-7(^P^z5[/(^](x; x') + iS[f^]{x- x')^(^7^) 



4i7^ rri-^i 



r(f+z^) 



(47r)T D 
(47r)f ^ 2 



r(i+.f) 



(92) 
(93) 



+ 



H ^ 
H 



^{^yM^+'^)+H^-'^)]+0(.)\g,,. (94) 



Because the differentiated fields in T^^ cannot contribute leading order log- 
arithms we see that the stress tensor takes the form, 



T, 



flU 



(95) 



where the potential is, 
H"^ fri 



-7 + C(3) 



, / + -(—) -y(-l)" C(2n-1)-C(2n+1 



(96) 
(97) 
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This is not the same potential V{(f) we found in section 3! Unhke that 
potential, Vs{ip) has positive curvature at ip = 0. However, the leading asymp- 
totic behavior for large cp is the same. 

The evolution of ip is controlled by V{(f). So the scalar rolls away from (p = 0, 
even though this initially means moving up the potential Vs{ip). 

What Vs{ip) gives is the expectation value of the operator which is the 
source of gravitational back-reaction. The two potentials disagree because 
the scalar effective potential V{(p) involves the square of the Hubble parame- 
ter, which is really — ^ R/12 for a general metric]^ That has consequences 
for the way V{(f) sources gravity. One can see this from the familiar case of 
a conformal coupling term in the Lagrangian, —^{p'^R^/—g. In 4-dimensional 
de Sitter the Ricci scalar is i? = 12H'^, so the contribution to the scalar 
effective potential is, 

AC = -l(^2^y3^ =^ AV{^) = GH'^if'^ . (99) 

But the stress tensor takes account of the way the Ricci scalar depends upon 
the metric for a general geometry, 

(100) 

Differentiated scalars cannot contribute leading order infrared logarithms, 
and the de Sitter Einstein tensor is R^^ — ^g^^R=—3H'^gf^^, so the induced 
potential in the stress tensor is, 

AC = -^(/?2^v^ =^ AFs(v') = 3^rV^ • (101) 

^This nonlinear dependence upon the Ricci scalar does not imply the usual kinetic 
instability associated with higher derivative gravity |47j . 
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6 Discussion 



Any theory which includes either massless minimally coupled scalars or gravi- 
tons will show infrared logarithms in the expectations values of certain op- 
erators. These enhance loop effects by powers of ln[a(t)], where a(t) ~ e^*^ is 
the inflationary scale factor. Of course loop effects are still down by powers of 
the (presumed small) loop counting parameter, but continued evolution must 
eventually bring about a situation in which the factors of ln[a(t)] overcome 
the small loop counting parameter and cause the breakdown of perturbation 
theory. Starobinskii has long advocated gaining control over this nonper- 
turbative regime by studying the series comprised of the leading infrared 
logarithm at each loop order [20] • In fact he has completely solved for the 
leading log limit of a massless, minimally coupled scalar with arbitrary po- 
tential [21]. 

In scalar potential models every field is "active." That is, it can poten- 
tially contribute to an infrared logarithm. Because the leading log result 
requires that all fields at a given order contribute infrared logarithms, one 
can perform an infrared truncation on the fields. This is why Starobinskii's 
formalism ends up being so wonderfully simple. More general theories also 
possess "passive" fields which cannot themselves contribute to an infrared 
logarithm. However, these passive fields can propagate infrared logarithms 
obtained from interaction with active fields. They can also mediate interac- 
tions between active fields. Differentiated active fields play much the same 
role. 

In propagating infrared logarithms, and mediating interactions between 
active fields, the ultraviolet parts of the passive fields (and differentiated ac- 
tives) contribute on an equal footing with the infrared. It is therefore invalid 
to infrared truncate either passive fields or differentiated active fields which 
appear in interactions. The correct procedure instead is to formally integrate 
out the passive fields, and the differentiated active fields, both from the action 
and from whatever operator is being studied. The expression which results is 
generally not local, but it contains only active fields. Because the nonlocal- 
ity is confined to inverse differential operators of passive fields, which cannot 
cause infrared logarithms, the associated Green's functions are always dom- 
inated by positive powers of the scale factor whose explosive growth weights 
the result completely at the upper limit in the leading logarithm approxima- 
tion. Hence the nonlocal effective action degenerates to a completely local 
and computable effective potential. At this point one has a local potential 
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model of the form Starobinskii has aheady solved in toto [201 1^ - 

Massless Yukawa theory is a wonderfully simple testing ground for these 
ideas because it contains a passive field — the fermion — without any dif- 
ferentiated active fields]^ It also lacks the subtle gauge fixing problems of 
SQED [25] and gravity [261 [IB]- In this paper we have exploited the classic 
solution of Candelas and Raine [21] to derive the Yukawa stochastic effective 
potential V{ip) (jSHD- We have checked the technique with an explicit two 
loop computation of the coincident vertex function. The result is in perfect 
agreement with the stochastic prediction. 

We have also obtained a leading log result for the stress tensor as a 
function of the scalar. Although this stress tensor takes the form —g^yVs{'^), 
our result for Vs{ip) (jHSD it is not quite the same as the effective potential 
V{ip) fl59l) which governs the scalar's evolution. The reason for this is that 
both potentials depend upon the dimensionless quantity, {fip/Hy, and the 
factors of = R/12 in this exert a nontrivial influence upon the way in 
which this model sources gravity. So one determines the scalar's evolution 
using V{ip), and one finds its impact upon gravity from Vs{(p). The two 
potentials differ, but they are each correct. 

A curious and potentially significant feature of both potentials is that 
they are unbounded below. The physics behind this seems to be very solid: 
inflationary particle production drives the scalar away from zero, which in- 
duces a fermion mass. That increases the magnitude of the fermion 0-point 
energy, which makes for a negative effective potential because fermion vac- 
uum energy is negative. We note that scalars seem always to induce growing 
mass [m [121 SHI US], so we expect that the effective potential of SQED will 
be positive for large fields. By comparison, gravitons seem to induce a grow- 
ing field strength renormalization [19j. It is intriguing to speculate on what 
that might mean for back-reaction in theories of gravity plus matter. 

Because the scalar effective potential is unbounded below, this model 
should decay forever. However, Hubble friction will make the evolution dom- 
inated by inflationary particle production until the scalar reaches nonpertur- 
batively large values. Although the initial effect is to raise the gravitating 
energy density, the large field results for V{ip) and V^(v?) agree. It seems 
inevitable that pure Yukawa must be unstable against slow decay to anti-de 

^Yukawa theory is so much simpler than SQED that this paper was complete well 
before a very similar analysis of SQED which was begun at approximately the same time 
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Sitter. In the Standard Model of flat space this same tendency is controlled 
by the positive effective potential from the gauge bosons. One naturally won- 
ders what the result might be for inflation, and whether or not this might 
be parleyed into a model in which inflation might be gotten to end without 
endless decay into anti-de Sitter. 



7 Appendix: Integrals from Section 4 

In section 4 we reduced the perturbative expression for the two loop con- 
tribution to ([n2D to a sum of dimensionally regulated integrals (|71|) . (17^ . 
( TTS]) . (|5U]) . (1521) and The next step is to partially integrate the inverse 

powers of Ax^ until they become integrable in D = 4. There is no distinction 
between ++ and H — terms at this stage. The identities we need are, 

1 _ ^2 / 1 N 

Ax3^-8 ~ {3D - 10)2(D -4) V Ax^^-io ^ ' 

1 / 1 X , , 

(103) 



Ax3^-6 (3D-8)(3L'-10)4(D-3)(D-4) VAx^^-iO' 

1 5^ / 1 , , , 

(104) 



Ax3^~4 (3D_6)(3D-8)(3L>-10)8(D-2)(L)-3)(L>-4) VAx^^-io. 

1 _ , 1 ^ 

A,t2^-4 ~ (2D - 6) (D -4) V ) ' (10^) 

1 9^ 



Ax2^-2 (2D-4)(2D-6)(D-2)(D-4) VAx2^- 

Because we are integrating over x'^, derivatives with respect to can 
be taken outside the integral, leaving an integrand which is flnite in D = 4 
dimensions. The limit D = 4 could be taken at this point except for the 
factors of 1/ (D — 4) which were picked up from the last partial integration. 
To segregate the divergence on a local term we add zero in the form, 

^2 ( - _ ^ = 92 f . (107) 

VAx:f-2^ r(f-l) ^ ' ^Ax^z^^ ^ ' 

Once this has been added we can take the limit D = 4 in the nonlocal term. 
We will work this out for the H — term, 

\ I \ _ \ 1 ^ 

2D-8 j Ax3^-io J ~ 2D-8jAx3^-io Ax:^r2 
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(108) 



3(D-4) 



2D-8 Aa;2 



^3(D-4)^2^ 



kD-4) 



-i(D-4) 



^ln2(^2^a;2_)+0((D-4)- 



Ax2 



(109) 
,(110) 



The other result H — we need is, 



f i_ 

L)-4 \ Ax2^ 



92 



/i 



D-4 



D-4 Ax2^-6 



2(D-4)^2^ 



-| ln(/i2Aa;2_) +^ ln2(/i2Aa;2_) +o((D-4)' 

A^2 



;iii) 
;ii2) 



It is only on account of the explicit factors of 1/{D—A) in Bi{x) and cot(|I?) 
in B2{x) that we must keep the order (D — A) terms in relations (11 101) and 
(11121) . The analogous ++ relations are, 



j 1_ 

2L)-8 \Ax3^- 
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2D-8r(f-i) 



ilnf^/i^Ax^ 



(Z>-4) 



ln2(/i2Ax2^)+0((D-4)= 



Ax"' 



..(113) 



D-4\Aa;2^-6j D-4r(f-l) ^ 
, ^2..-4),2 f -I + ^ M VAx^,,) +0 ((D-4)^, ^^^^^ 

The use of these partial integration identities results in each of the terms 
dZD, dZHD, dZH]), ([80]), ([82]) and dHU) having a finite, nonlocal part and a 
potentially divergent, local part. For A{x) these are. 



,-ffl-4 . rlnfu^Axi^) Infu^A^ 



' irO (3£>-6)(3i5-8)(3£>-10)(2D-4)(2B-6)(2D-8)'^ ' 
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Before giving the nonlocal and local terms for the Bi{x) and Ci{x) we will 
evaluate the two above. The local term is simple because the derivatives in 
the numerator are so easy, 



2 9 a 



(2-f)(3-f)(44)(5-f)^V 



8-213 



(117) 



The factor of (2 — means that Al{x) is actually finite and we can set 

D = 4, 

M^) ■ (118) 

To evaluate any of the nonlocal terms it is best to extract two more 
derivatives, 



2i337r6 



ln^(/i2Aa;2^) -2ln(^I^^Ax 
ln'(//2Aa;2_) +21n(/x2A 



2 

++ 

2 

X 



(119) 



The point of doing this is that we can now exploit the exact cancellation 
between ++ and H — terms outside the past hght-cone. (Note that we do 
not want to do this before the last derivative is extracted because the limits 
of integration must be constant for us to extract derivatives.) We define the 
temporal and spatial intervals in the obvious ways. 



A77 = 77 — 77' 
It is then straightforward to show. 



Ar = \\x — x' 



(120) 



Ini^fi^Axl^j - \n[i2'Axl_j = 27rie{Ar] - Ar) , (121) 
In2(^2^x2^) - ln2(^2^4_) = 47ri^(A77 - Ar)ln[fj,\Ari^ - Ar^)]. (122) 
These relations bring the nonlocal term to the form, 

An{x) = ^n^^^ fd^' J d'^x'e{A7]-Ar)lln[i/{Ar]^-Ar' 



1 , (123) 



where the initial time is 77/ = —1/H. Note that Aisf{x) is now manifestly 
real. 

The next step is to make the change of variables f—x' — x and perform 
the angular integrals, 

-fa-' 



An{x) 



f)8 

2937r4 



d7]' 



VI 



(irr^ I In 
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lj:'{A7i'-r') 



- 1 



(124) 



We then make the change of variables r = Arjz and perform the integration 
over z. 



An(x) 



-fa-' 
2937r4 



9 



(125) 



Owing to the factor of Ar/'^, three of the external derivatives can be brought 
inside the integral, 



-fa-' 
2937r4 



rdr]'A\n{2fiAr]y 



(126) 



At this stage one performs the integral over r]' and acts the derivatives, 



AM 



fa 



2737r4 



dUAT^, ln(2/xAr/,) - 1 



fH' 



2H^ {a-iy 



(127) 



For the 5-terms it is best to convert from D to e = A — D. All three of 
the B terms (ff6l) . (j78ll and (jSOll . contain an overall factor of. 



2^7rt^a^-i 



^ r(2-|)r(3-e) 

24-^7r6-f^a3-^ 

The integrand of the Bi term fl76l) is this overall factor times. 



(128) 



/I 2 J 1 

-H aa a < 5- 



D-4:T{D-1) 

[\H^aa'Y'a'd^ 
^(4-3e)(2-3e)(2-2e)(-2e) 



1 



3D- 

++ 



Ax 



4-6^ r^(2-f) 
r(3-6) 



eln(/i2Ax^^) + ^^(/i^Ax^^) + . . . 



+ ^^^^Y^5^(^-^')^ (129) 



The associated local term is, 

fH^^^-^^ (l-|)(l-|)r2(2-f) 



327r4- (l-|6)(l-e)(l-|6)e2 

/gV"^-(l-fi)(l-f)r^(2-f) 
16vr4- (l_3e)(i_e)(l-|e)e2 



,2e 



(130) 



(131) 
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The integrand of the B2 term fl78p is the factor 01281) times, 
1 1 



7r 



— vr cot(^—D]a' 



/ 71 e \ 

TTCOtf^ — j 



X 



(4-2e)(2-2e)(2-e)(-e) 



ln(/i^Ax^J + |e^ln"(/x^Ax^^ 



The appropriate local term for B2 is 

B2l{x) = 



/i72-.^-7rcot(f)r(3-6)^_3^,^2^ 



4-e 



cot 



2 

7re 



The integrand of the term (|80l) is the factor (11281) times, 



a' In(aa') | 



1 



a' In(aa') d"^ 



2 A ^2 



ln(/i^Ax;^) + . 



(4-2e)(2-2e)(2-6)(-6) 

^2-ie,,-e 



The appropriate local term for B^ is. 



(132) 



(133) 
(134) 



B3l{x) 



3-^7r4- (i-|)(i-e)e 



{ln(a)92a + 92(aln(a))} , (136) 



/ff4-/i-r(l-e) 



57r 



4-e 



{ln(a) + -}a' 



(137) 



For the nonlocal terms it is useful to extract a factor to go with ([12 
making a total multiplicative factor of. 



_,^^.-r(2-|)r(3-e 



2'-'n'-i'a'-' 16(l-|)'(l-e) 



27-^7r6-^a3-^ 



(138) 
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It is also useful to note the expansions for the Gamma function, 



r(i-^) = 1+76+ 

and the cotangent, 



7" tt" 
¥"^12 



r^(i-f) 
r(i-e) 



TT 



(139) 



TT 



cot 



1 



— ( 

12 



0(e 



(140) 



We can expand the nonlocal integrands without regard to ± variations. 
The integrand for the nonlocal part of -Bi(x) is (11381) times. 



(l-t)(l-ff 



r2(i- 



(l-|e)(l-e)(l-|e) r(l-e) V 



1^ 
e 



^,^ H^aa\ ^d^j ln(/i^ Ax') + e In' (/i^ Ax') 



ln(/i'Ax" 
Ax2 



ln(;u'Ax' 



Ax2 



Ax2 

rin'(/i2Ax') 



ln(/i'Ax' 
Ax2 



Ax2 



-0(e). (141) 



The integrand for the nonlocal part of B2{x) is fll38p times, 
TT /vrex 4 1 ln(/i'Ax') + |eln'(/i2Ax') 1 



Ax2 



Ax2 



4 



ln"(/i2Ax^ 
Ax2 



+ 0(6). (142) 



And the integrand for the nonlocal part of B^{x) is f llSSp times. 



ln(/i2Ax' 
Ax2 



+ 0(e). 



(143) 



Combining (1141111431) . multiplying by (I138p . including the integral, and taking 
e = 0, gives the following result for the nonlocal part of the B term. 



Bn{x) 



'2^ 



a-^d^ d^x'a' 



+ ln(|)l ln(^2^x2j + iln2(/.2Aa; 



Axl, 



ln(/i2Ax^_) + |ln2(/i2Ax^ 



• (144) 
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The reduction of Bjq{x) proceeds very much hke that of Aj^{x). We 
begin by extracting another d'Alembertian, then combine ++ and H — parts 
to make causahty and reahty manifest, and perform the angular integrations, 



Bn{x) 



a-^d^ d^x'a' 



^ln^(/x2Ax^J + |ln(f ) ln^(/.2Ax^,J-iln(f ) H^^^^xl 
[-^ln^(/.2Axt)-|ln(f)ln^(/.2Axt) + iln(f)ln(/.2Ax^ 

iln^[/.^(Ar?^-r^)] -g+ln(^) ln[;.^(A^^-r^)] -ln( 



, (145) 



2^ 



(146) 



The next step is to make the change of variables r = ^rjz and perform the 
integration over z. For this purpose it is useful to note the integrals, 

j^dzzH.{-^)=-- and y^d..Mn^(^) = ---. (147) 
Applying these identities gives. 



17)1 1 3 

--ln(2;.A,)+--- + ln(-)[-ln(2;.A, 
The next step is to bring three of the derivatives inside, 
Bn{x) = -^^,a-'dlj\'a'^rf{\n\2i,^ri) 



11 



;i48) 



21n(2/xAr/)+2-y+ln(^) [21n(2/xAr;) - 3 



2// 



(149) 



2^4 ^ 



a'^d^J\' a' At] {2 ln2(2/iA?]) 

-21n(2//A77)+2-^+ln(|-) [4\n(2fiAr] 



(150) 



—d^J^dr]'a' I 21n2(2/xAr;) +41n(^) ln(2/xAr/) -y i . (151) 



2%^ 
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Before performing the integral it is best to rearrange the integrand, 



3 „ 2 



ln^(2/xA7?)+21n(^) Inh/iAr))-^ = ln^{HAr))+3ln{HAr))+K , (152) 
where the constant K is, 

i..-.„^(|)+31n(|)-f (153) 

The integral gives a complicated result, 

i?iv(x) = ^a-=^ao^|--ln-^(a)--ln(a) + 2C(3)-25: — 
-l„(a).n^(l-i)+2.,>(l-i)i:(^ 

n=l ^ ^ „=i ll- ) 

In acting the derivatives there is no point to keeping any but the logarithmi- 
cally enhanced terms, 

Bn{x) = Jl-a-'dl[-l\n\a) + ^\n\a) + 0{ln{a))Y (155) 

= -^K(«) + 0(1)}. (156) 

By comparison the C terms are straightforward. The two nonlocal terms 
cancel. 

The corresponding local terms are, 

c (X) - I /-^V-" (i-i)r(2-i)r(4-|) a^^ 

- ^^^^ (l-|6)(l+6) V ^^^^^ 

^^^^^^ - "26-7r^- (l-t)(l-6) 7- 
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